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We examine the nuclear liquid-gas phase transition at large number of colors (Nc) within the 
framework of the Van Der Waals (VdW) model. We argue that the VdW equation is appropriate at 
describing inter-nucleon forces , and discuss how each parameter scales with Nc. We demonstrate 
that Nc — 3 { our world ) is not large with respect to the other dimensionless scale relevant to 
baryonic matter, the number of neighbors in a dense system A^'jv. Consequently, we show that the 
hquid-gas phase transition looks dramatically different at A^c — > cxd with respect of our world: The 
critical point temperature becomes of the order of Aqcd rather than below it. The critical point 
density becomes of the order of the baryonic density, rather than an order of magnitude below 
it. These are precisely the characteristics usually associated with the "Quarkyonic phase". We 
therefore conjecture that quarkyonic matter is simply the large Nc limit of the nuclear liquid, and 
the interplay between Nc and Nn is the reason why the nuclear liquid in our world is so different 
from quarkyonic matter. We conclude by suggesting ways our conjecture can be tested in future 
lattice measurements. 
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I. INTRODUCTION 



Strongly interacting matter at moderate (~ the confinement scale) quark chemical potential and moderate 
temperature T has recently received a considerable amount of both theoretical and experimental interest. Such 
matter can hopefully be produced in heavy ion collisions and is thought to exhibit a rich phenomenology. The 

latter includes one 5] or more @ critical points, spinodal instabilities precursors to color superconductivity B , 
separation between chiral symmetry and confinement chirally inhomogeneus phases [isl - flq . new phases [Toj 

etc. 

These conjectures are, however, extraordinarily difficult to quantitatively explore in a rigorous manner. The quark 
chemical potential fiq is nowhere near the asymptotic freedom limit where perturbative QCD can be used [2^. It is, 
however, way too high for existing lattice-based approaches, dependent on fiq/T <^ 1, to work [2l| . 

Perhaps the only relevant quantity with can be uncontroversially be called "a small parameter" (albeit not so small 
in the real world!) is 1/Nc, where Nc is the number of colors '22']. While the large asymptotically Nc theory shares with 
QCD asymptotic freedom for hard processes and confinement for soft ones (separated by the energy Aqcd ~ 250 
MeV, independent of Nc) , the Nc scaling of different observables can be used to establish a model-independent 
hierarchy. Thus, the shape of the phase diagram can be said with relative certainty to look like Fig. [TJ Phases I 
and III are, respectively, the familiar confined chirally broken Hadron gas (where pressure ~ iV") and the deconfined 
chirally- restored quark-gluon plasma (where pressure ~ N^). Since at large Nc gluon loops dominate over quark 
loops, the critical □ temperature Tc ~ iV^Agcc and the critical chemical potential necessary for deconfinement is 
very high, fig N^Aqcd- 

Consequently, in the large Nc limit, the phase transition line becomes horizontal for moderate iiq. In this limit 
the transition between zero baryonic density and finite baryonic density matter is infinitely sharp at NcfJ^q ~ 
NcAqcd ,19], since the baryon density ^ exp[—Nc{AQCD — Mg)] goes to zero exponentially with Nc for chemical 
potentials less than the baryonic mass. Thus, a new phase (II) emerges where the nuclear density is 0(1) Aq(^^, 
parametrically much less then that required for deconfinement, O (Nc) Aq^j^. 

Naively, since fiq ~ Aqcd is nowhere near the chemical potential required for deconfinement, this phase should just 
be that of dense nuclear liquid (the large Nc limit of the nuclear liquid, well-studied theoretically and experimentally 
(23 - [30| ). where nucleons are close to touching each other, yet confinement is still there and degrees of freedom are 
baryons and mesons. The critical parameters for this transition, however, are far above what is seen in the real world, 
in line for the much stronger nuclear force seen at large Nc wrt Nc = 3 [3l|. Moreover, as pointed out in T^, at 
this chemical potential inter-quark distance ~ ^/Nc, leading to the apparently paradoxical situation of quarks close 
en oug h to interact perturbatively in a confined medium. 

[l9| proposed to solve this conundrum by postulating that in the new phase the quarks below the Fermi surface 
act as free objects but the Fermi surface excitations are confined. Thus, while the new phase is confined, the entropy 
density and pressure feels the quark degrees of freedom and ~ N^ rather than ~ as in the usual hadron gas. This 
new state of matter, called quarkyonic in |19l] , should also be realized in our Nc = 3 world and reachable in heavy ion 
collisions since large Nc is at least qualitatively true in our worlcfl. 

A great deal of investigation has gone on to see weather quarkyonic matter appears in any effective theory of QCD. 
While a phase transition does seem to exist which has some of the characteristics described above ^12hl4j , it is not 
clear weather the most interesting properties (P ~ Nc and chiral symmetry restoration in the confined medium) are 
physically realized, as we do not have a model realistic enough but still computable. Other approaches have found 
no evidence for a ny such transition [33. , .34J ,or have claimed the "quarkyonic" phase to have different properties for 
those claimed in [l9| (eg [35[ conjectures a chirally broken but deconfined constituent quark plasma). 

As discussed in the introduction, the main difficulty of theoretical investigation in this regime is that there is no 
reliable approximation technique which is capable of distinguishing between models. We are thus left with effective 
theories, such as the NJL and pNJL model [l3j [13, [HI , or Gribov-Zwanziger confining gluon dynamics [lB|- The 
results obtained with these models, however, are highly dependent on the assumptions made in them, assumptions 
which can not be rigorously shown to derive uniquely from QCD. In case of the critical point [36[, different models were 
shown to give very different answers. Additionally, none of these models contain features unique to non-perturbative 
QCD, such as exact quark confinement. As a consequence, the crucial aspect of the quarkyonic hypothesis, scaling of 
entropy density with Nc in the quarkyonic phase, can not be adequately tested with models such as pNJL (T^. [l3|. 

A possible way out are techniques deriving from Gauge-string duality [37| . While no string theory with a dual 
looking like QCD is known, several models were developed which share with QCD some of its more notable non- 



^ Unfortunately the standard nomenclature is somewhat confusing as the subscript "c" means "colors" in Nc and "critical" for thermo- 
dynamic quantities (the critical temperature Tc, density pc, pressure Pc, chemical potential fic and so on) 

^ Recent work has broadened the definition o f "q uarkyonic matter" to other characteristics, in particular related to phenomenology of 
chiral symmetry breaking and restoration Illl4l6ll . In this work we use "quarkyonic matter" to refer to matter which is confined but 
whose pressure scales with Nc,in accordance to the definition given in [19| 
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FIG. 1. (color online)The phase diagram for large N^. See text for a description of the phases I,II,III in various models. See 
footnote [T] about the nomenclature 



perturbative characteristics, such as confinement and chiral symmetry breaking 38]. These models can be used to 
extrapolate to regions inaccessible to pQCD and the lattice, while retaining qualitative aspects of non-perturbative 
QCD such as its strongly coupled nature and dynamical confinement. 

A finite chemical potential study l39l. liol within the Sakai-Sugimoto model [s^ has shown that the system looks 
similar, but is crucially different from [l9| in several ways. The basic structure of Fig [1] which emerges is common 
to [l^: The temperature between phases I (confined) and III (deconfined) is insensitive to N^, ~ -/V°. The 
curvature of the deconfinement transition line wrt the chemical potential becomes negligible in the large Nc limit 
(The transition to QGP in the jis plane grows as A^^ )• Most intriguingly, a new phase II emerges, with the transition 
line at /ig ~ O (1) N^Aqcd and the nuclear density as the order parameter, just like in [T^. 

There are, however, profound differences: (39| finds that both phases I and II are confining and chiral-broken. No 
evidence exists that the scaling of the pressure changes between I and II. In fact, the only difference between I and II 
seems to be a discontinuity in the Baryonic density. The authors of ,3.S] interpret phase II as the well studied nuclear 
gas liquid phase transition (23l - [30| , rather than as a new undiscovered phase. If this interpretation is correct, than 
searching for the quarkyonic phase and/or the triple point separating I, II, III at upcoming low energy experiments 
[ll-ll] would be fruitless, as in our Nr = 3 world the liquid-gas phase has been extensively studied theoretically j23l - [27{ 
and pinpointed experimentally |28l - l30| . and its transition line is understood to lie well below Tc, so that no triple 
point exists. 

This ambiguity on the phase diagram is compounded by a limitation of our understanding of baryons and their 
interactions in the large Nc limit. Baryons at large Nc have long been thought to be well-approximated [3l| by the 
solitonic "Hedgehog" configuration. Fig. [21 A "flower" of Nc quarks linked to a Baryon junction via strings (the 
confining potential). The bulk of the baryon mass is carried by the confining strings, which therefore ^ Nc, with 
the proportionality constant given by the quark mass or the string tension. The baryon radius, on the other hand. 

It is generally believed that in the large Nc limit baryons, unlike mesons [l^ are not weakly bound states [31], but 
that 2-baryon, 3-baryon and A^-baryon forces generally ^ Nc- The comparatively weak nuclear force in our world is 
therefore the result of an accidental cancellation of attractive and repulsive forces. 

Such a scaling, however is somewhat counter-intuitive given the natural hierarchy between multi-body forces arising 
in effective field theories [4l|,|42]. Furthermore, it was found [i^ that the Nc dependence becomes faster at higher 
order, [i^ casting doubt on the existence of a well-defined large Nc hmit of nuclear forces. Recently, the picture has 
become even more controversial with the conjecture [44| that the consensus of the skyrmion picture of baryons at 
large Nc is flawed, and that consequently the nuclear force scales as A^°. The rather weak inter- nuclear potential (in 
comparison to the forces holding quarks and gluons together) observed at A^c = 3 in this picture would therefore be 
natural rather than accidental. 

This theoretical ambiguity leaves space for exploration with simple effective established models. Perhaps the oldest 
and most well-known model capable of capturing the high-temperature dynamics of the non-relativistic liquid-gas 
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FIG. 2. (color online)A gas of baryons in the large Nc limit interacting 



phase transition is the Van Der Waals (VdW) equation of state in terms of Pressure P, Temperature T and nuclear 
density p . This model is undoubtedly extremely rough, and can not be trusted beyond the level of a qualitative 
estimate. It does, however, have the virtue of being universal, a mean field expansion general enough to describe 
systems as different as water-vapor, nuclear matter[45i| (where the potential has been successfully related to effective 
theories of QCD and charged black holes 0,111], the latter generically being the dual of strongly coupled thermal 
Gauge theories at finite temperature and chemical potential. 

Changing the form of the Baryon-Baryon interactions will not change the form of the Van-Der- Waals terms (or the 
higher order Virial corrections), but merely change the numerical values of these terms. Hence, in the vicinity of the 
phase transition, this model can be trusted to provide a good qualitative description of the relevant physics (such as 
scaling with Nc), though corrections of O (1) at any given Nc, to be calculated by more sophisticated models or on 
the lattice, would be expected. 

We also note that the VdW model has qualitative limitations which prevent it from modeling the "right-hand side" 
of the phase diagram (the dense liquid phase), to check, for example, if in it the pressure would ^ Nc (as proposed 
in flo'l for the quarkyonic phase) or ~ N^ (as is in our world for the nuclear liquid). This is because in this model 
interactions are integrated out to a modification of the dispersion relation of nuclei, rather than treated as additional 
degrees of freedom that can also contribute to pressure. This is why the VdW equation of state fails to describe the 
electron pressure in a crystal, and why it would be inappropriate to describe the quarkyonic phase, where pressure 
is dominated by partons trapped deep inside the Fermi surface [l9{ . These limitations, however, do not prevent a 
qualitative estimate of the phase transition line of the nuclear gas, they just prevent a peek into what happens after 
the system crosses that line. In this work, we shall content ourselves with the first task. 



II. THE VAN DER WAALS NUCLEAR GAS AT LARGE Nc 

In the large Nc limit, the only A^'c-invariant scale of the theory is Aqcd, the scale at which the 't Hooft coupling 
constant becomes A ^ O (I). While a precise value of this scale depends on the scheme used to calculate it, its roughly 
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Aqcd ~ - 200 - 300 MeV That this value is very close to the "constituent" quark mass, ms/Nc ~ 310 

MeV, is not so surprising given the relationship between confinement and chiral symmetry breaking [50| : Assuming 
that confining effects set in at the scale where A > 1, and solving for a system of massless fermions in a confining 
potential (such as a simple harmonic oscillator) will give a "dynamical mass term" (or, equivalently, a breaking of 
chiral symmetry) of the correct order. Physically, the constituent quark masses will be given by the inverse of the 
baryon size, which is set by Ag^^. 

Since the only scale in our theory is AQC£i,it is natural to expect that any physical quantity is ~ /(7Vc)Aqp^, 
a dimensionless function of Nc times a power of Aqcd set by the dimensionality d of the quantity. Henceforward 
we shall adopt this assumption, and, for brevity, set Aqcd to unity in the equations. The reader should multiply 
any dimensionful quantity in the equations by the appropriate power of Aqcd (For example, the Baryon mass is 

NcAqcd in the text, and ~ Nc in the equations). 

In this notation, the Van Der Waals parameters a, b and the curvature correction become dimensionless a,/3,7 times 
the appropriate power of Aqcd (3 for a, 2 for /3,4 for 7), and the VdW equation Q becomes 

{p-'-a){P + pp^-jp^)=T (1) 

Here, a parametrizes the excluded volume and /3,j the interaction. When a,/3,^ = the system reduces to the 
classical ideal non-relativistic gas. In the VdW equation usually referred to in textbooks, 7 = 0. This system is 
solvable more easily, but, as can be noted, at T = P is generally less than zero. The next-to-leading order, 7 > 
"curvature correction" fixes this artifact, and becomes dominant at low temperatures and high densities. 

The VdW (7 = 0) equation can be shown to arise as an effective first order description of any interacting non- 
relativistic system where the interacting potential depends on 2-particle interactions only, and the system is dilute, 
so that the minimum distance between particles is far above the excluded volume. Note that a 7 term arises as a 
correction to either of these assumptions, since both 3-body forces and a 2nd term in the Virial expansion can be 
shown to give rise to 7 > 0. Higher-order corrections behave in a similar way. 

The first assumption of the VdW set-up (2-body interactions dominate) seems to be true in effective field theories [4ll 
[in, since the strength of an n— body interaction typically ~ (/c/A)" where k is the relevant momentum scale (e.g. the 
binding energy) and A 3> A: the scale at which the effective theory breaks down. 

We note that this reasoning is not in contradiction to (sij , since its perfectly possible for n-baryon interactions to 
~ Nc{k/A)"-. Thus, 2-body forces are more important than three body forces which are more important than 4-body 
forces and so on, but all scale equally with Nc- From this, we learn that generally 7 < /3 but their Nc dependence 
should be the same. 

The second assumption of the VdW set-up (inter-nuclear distance ^ excluded volume) is certainly true in our 
world (the inter- nuclear separation is roughly an order of magnitude larger than the nuclear size), but not necessarily 
in the large Nc world. To understand why this is the case, we can see Fig. [5] (in the spirit of the nucleon hedgehog 
model presented in [3l|): In our world, packing hadrons in such a way that there is a significant overlap area would 
result in a deconfinement phase transition; The percolation picture of deconfinement [sil . [s^ makes this physically 
intuitive. 

In the large Nc world, however, the fraction of quarks which overlap becomes negligible (~ Nc), while the bulk 
of the baryon mass remains in the baryon junction. It is therefore possible for baryons to start overlapping (nuclear 
density reaches ~ Ag^^^) while remaining well-defined objects. 

As the number of colors — 00, it is therefore logical to postulate a will reach the limit of a ^ O (1) Ag^^. At 
Nc decreases from cx3 to O (1), a can only increase, so a useful ansatz, universal in any limit where a can be Taylor- 
expanded, is a ^ (1 + A/Nc) Aqcd- To understand the behavior of A, it must be remembered that it must sensitive to 
N]\[, the number of neighbors a nucleon has in a tightly packed nuclear material. The more neighbors, the more Pauli 
blocking of valence quarks must be important, and the more the presence of neighbors will disturb the configuration 
space part of the quark wavefunction inside the nucleons. Since, due to confinement, any such disturbance of the 
nuclear wavefunction adds an energy of ^ Aqcd, the nuclear repulsive core will be larger than the inverse of the 
nuclear separation up to the deconfinement temperature. If the number of colors is larger than Nn, this problem will 
not exist since it will be possible to arrange quarks so neighbors will be of different colors. In this limit, therefore, 
baryons can be tightly packed (interbaryonic separation ^ Aqcd) without the configuration space part of the baryonic 
wavefunction being disturbed. 

Nn, of course, is a function not of Nc, but the (fixed) number of dimensions d and "packing scheme" , Nn ^ k{d)N^. 
The exact form of the "kissing number" function k(d) in arbitrary dimensions is unknown [53|, but seems to be 
approximated by k{d) ^ 2°"^, with a transcendental a — 0.22... . /c(l,2,3) is, respectively, 2,6 and 12. 

Therefore, we assume that 

"-(f) — — 



(2) 
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In a sense, the relation above is a parametrization of experimental data combined with natural constraints: We know 
that the nuclear size can only go as ^ A + B/Nc. We know that in our world the inter- nuclear distance is ^ ^qcd, 
but it can only decrease with Nc. Putting the ~ term to O (1) and the ~ term to O (10) is the only way to 
account for these limits, as the results of the next section shall show. The interpretation due to the kissing number 
is then an appealing physical explanation. 

Note that in our world the first term dominates, since for any realistic packing of nuclear matter the number of 
nearest neighbors is considerably larger than Nc — 3 (eg for cubic packing, its 6 or 8, depending on weather corners 
are included. Generally it should go as k{d) — 9,parametrically larger than 3). Pauli blocking between valence quarks 
of neighboring nucleons, therefore, can not be ignored and keeps the inter-nuclear distance parametrically larger than 
the nucleon size. In the large Nc but three-dimensional world the second term takes over as the effect of Pauli Blocking 
becomes negligible, and the excluded volume approaches Ag^^, since baryons become more tightly packed in the large 
Nc limit (presumably, the "jamming" phase transition described by percolation [5l|, [s^] would be a good description 
of the liquid-gas transition in this regime) . Note that the small parameter behind the VdW expansion (the scale of 
the excluded volume over the scale of the inter- nuclear potential) is < 1 irrespective of Nc-, even through it reaches 
1 asymptotically with Nc. For a rough qualitative estimate, therefore, the VdW equation is always justified, but we 
have to keep in mind that higher order quantitative corrections become increasingly important as Nc — oo. These 
corrections are known to influence the shape of the phase transition line rather than its basic limits. 

The behavior of /3 (and hence 7) is somewhat ambiguous: As we remember from [54|, /3 is related to microscopic 
physics via 

/? = 27rr / drr"^ ( 1 - exp 
where V{r) obeys a class of Yukawa- type potentials 

Vir) ^ 'j^EtM (4) 
r 

converging to Coloumb as M — >■ 0. 

As pointed out in [sij, baryon-baryon interactions by both gluon exchange and meson exchange ~ Nc. This would 
mean that for a low nuclear potential {{V{x)) /T <C 1), /? ~ Nc It would seem likely, therefore, that the large Nc 
limit the nuclear "liquid phase" is actually a tightly bound ( solid?) material, where the VdW approximation stops 
holding. The behavior conjectured in [i^l, on the other hand, would mean that /? ~ Nc or ~ logA^c, allowing for a 
nuclear liquid phase not too dissimilar from our world. 

In this work, we shall test the consequences of both assumptions, by parametrizing /3 Nc, where v can be or 1 
{(3 ^ log A^c has the same qualitative dependence as Nc, through the convergence rate is parametrically slower). We 
shall assume that /3 and 7 have the same dependence on Nc, as expected in both iSlJ and 44i]. 



Vir) 
T 



(3) 



III. THE CRITICAL POINT 

From textbook formulae, [5^1 we can immediately read off the conditions (temperature Tc and density pc) of the 
critical point in case 7 = 0. As argued in the preceding paragraph, the critical point density 

should be <^ ^qcd '^^ world, but ^ J^qcd ^he large Nc world. 

Assuming /3 ~ Nc yields the following evolution for the critical point temperature 

Tc^ —^^[ — (6) 

27a \Nn^-Nc) ^ ' 

For the case of 7 ^ Nc ^ 0, the formulae become more complicated, but still analytically tractable, by solving for 

dP _ d^P 
dp dp- 

with the equation of state being an additional constraint. We obtain 

24iV,4 + 4jV,^iV^Fi + 2^NcNlD ~ m^F, -f SiV^f^ .,2 ^ 

288{Nc + NnY " ci/U cj K ] 



,2 - (7) 
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FIG. 3. (color online)The evolution of the critical point temperature (panel (a) ) and density (panel (b) ) with Nc in different 
scenarios, with 7 = (solid) and 7 ~ N" (dashed). Note that the density does not depend on v. Nn is assumed to be 10, the 
order of magnitude of the real world value. Temperature is presented in units of Aqcd, while baryon density in units of ~ 1 
baryon per fm~'^ 



where D = ^ ^Nl + 47VciVjv + 3iV^,i^i = (V3-D - STVat) and F2 = {-.fiD + GA^at) and 51,2 (^c) are rational functions 
of iVc where the powers of the numerator and the denominator arc the same. 

The evolution of the and pc are shown in Fig. [3] for A^at = 10. As we can see, the presence of v makes no 
qualitative difference for the critical point density pc, as generally believed and expected, pc goes to its asymptotic 
"overlapping nuclei" value regardless of v. The same, however, can not be said for Tc. 

For = we get the behavior corresponding, physically, to the scenario where the quarkyonic phase and the nuclear 
matter phase coincide: In the Nc ^ iV^r regime, ^ ^qcd-, while in the large N,. limit it saturates to ~ ^qcd as 
in [1^ and [39], with a numerical factor of 8/27 ior v — Q and ^ 0.31 for v = 1. . The numerical factor could very well 
rise to unity once a more realistic model (including excited nucleons or higher order corrections) is employed, through 
the deconflnement transition precludes it rising above unity. 

For V = 1, limAT^^^ooTc ^ NcAqcd- Since Aqcd ^ N^, this seems to signal that in this picture the the VdW 
scenario is inappropriate to describe the nuclear liquid-gas phase transition at large N^, as the critical point would 
soon overtake the deconflnement temperature already at ~ 10, and the binding between nucleons soon overtakes 
excluded volume effects (As conjectured in [l^ and revised in [11]). The applicability of VdW would therefore be an 
accident of our low Nc world, where the binding degrees of freedom do not contribute significantly to the pressure 
until well above the nuclear liquid-gas phase transition. 

Independently of the value of v, the location of the critical point in the real world, T <^ Aqcd,P ^ ^qcd^ can be 
understood by the realization that we do not live in a large Nc world, in the sense that Nc is not larger than any other 
dimensionless scale of the system: In our world the number of nearest neighbors dominates over Nc, lowering the 
scale at which the "dense" (liquid) phase appears to a density much lower to the overlap density of the hadrons. The 
evolution of the critical point from Nc ^ Nn to Nc ^ Nn tends to agree with the "nuclear matter" interpretation of 
the new phase as suggested by the calculations in ^ , though we can not rule out that this phase also has quarkyonic 
properties. 

In the next few sections we shall demonstrate that this conclusion is valid for the phase diagram shape, and not 
just for the critical point. 



IV. THE PHASE DIAGRAM FOR 7 = 



A. T — p plane 



We now examine the /? ~ iV^ case further, by calculating the density jump and the curvature in phase transition 
space. For simplicity, we shall concentrate on 7 = and leave the low temperature regime for a future work. As 
shown in the previous section, away from T ^ ps this approximation should be acceptable. 

We rewrite Eq. [T] as a cubic equation. 



a/3p3 - (3p^ + p{T + aP) ~ P = 



FIG. 4. (color online) The phase diagram on the T — p plane /3 ~ A^^ ( panel (a)) and P ^ Nc (panel (b)) for 3,5,8,20,100 colors. 
Higher curves correspond to increasing Nc. The dashed line in panel (a) corresponds to the asymptotic limit, absent in panel 
(b) (As in that case Tc ~ NcAqcd)- Nn is assumed to be 10, the order of magnitude of the real world value. Temperature is 
presented in units of Aqcd, baryon density in units of ~ 1 baryon/Jm"^ 



The task of finding the coexistence line of the phase transition diagram is a non-trivial one, giv ing rise to several 
computational approaches. In this work, we have used the parametric solution described in [55|. Given that the 
entropy difference between the liquid and gas phases is As, the solution can be parametrized by 



^±As/2 



'9,1 



a 



f 



As 



(10) 



where, pg i is the density in the gas and liquid phases respectively and 



fiy) 



y cosh y — sinh y 
sinh y cosh y — y 



The temperature can be found from the same parameters and the requirement that pressure in the liquid and gas 
phase transitions has to be the same, solving Eq. [T]for pressure, equalizing and doing some algebra gets us 



1 - 



Pi 



Pg - a 



(11) 



As can be seen in Fig. 2] (panel (a)), as long as the nuclear force potential V{r) ~ iV°, the VdW transition has 
a sensible limit. As Ai'c — )• cjo ,Tc — >■ O {1)Aqcd, Pi — Pg ^ BaryonxAgp^ Our world's parameters, where the 
critical points and the latent heat ^ Aqcd are also qualitatively reproduced, but they lie squarely in the "low Nc" 
{Nc < Nn) region. 

As expected, introducing a /? ~ Nc dependence makes the phase diagram's height shoot up to infinity, while keeping 
the width constant, an unphysical limit (Fig. Impanel (b)) 

Finally, we note that at large Nc the mixed phase occupies the bulk of the space between p — and p — Aqcd/3. 
Physically, this is reasonable, since in the diluted phase p ^ exp [—Nc] — >■ 0. Thus, adding "a few baryons" into 
the system while p <^ ^qcd ^^^^ ^'-'^ change the chemical potential. This indicates the mixed phase should go from 
p <C Aqcd to /9 ~ ^QCD^ exactly as our calculations show. 

Our world is, however, very far away from the large N- limit in this respect, since the nuclear liquid-gas mixed 
phase actually occupies a small area of the T — p plane [2J, [2^, [l^. As Fig. Fig. 2] demonstrates, the scaling of the 
mixed phase is well-accounted for by the interplay between the Nc ^ Nm and Nc ^ Nm limits. 



B. T ~ plane 



The chemical potential can be obtained [45l. [54. ISq. [57| by the textbook thermodynamic relation p — {dP/dp)j,. 
Inverting, and writing in terms of p,q = fiB /Nc we have 



Nc 



f{p',T)dp' + F{T) 



(12) 
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where the first term is the nucleon mass and 

Note that the integral has a logarithmic divergence at p ^ 0. This diverge is canceled out in Eq. [T^lby F{T), via 
the requirement that fiq{T, p — >• 0) — (So F{T) — —f{T, p — O),also a logarithmic divergence). The resulting is 
always well-defined, and equally valued for pg and pi. 

The ideal gas limit can be obtained by putting a = /3 = in Eq. [131 In this limit, f{p) = T/p and F{T) = 
InA"'^ where A is the thermal wavelength = ^J'rr^rJ^ ^ ^/N^T. We therefore recover the ideal gas formula, 
/iq = 1 + Tin (A^p). Note that, due to the scaling given by Eq. [2J a = is not realized at any Nc, since the excluded 
volume in a confining theory has to > Aq^^. In our world, however, one can generally neglect the excluded volume 
[srj as a at the liquid gas phase transition. This is however not true at Nc » Nn- 

At low temperatures a correction due to Fermi-Dirac statistics is necessary. Neglecting the effects of Fermi-Dirac 
statistics on the excluded volume and interaction, we would get that /i^ — > /i^ + AiipD, with the correction A/ifd is 
given by 

^^^FD = ^ [log(z) - log (pA^)] (14) 

, where the fugacity z can be obtained by solving the the equation linking the density of the ideal gas to the chemical 
potential [s^] 

Note that this correction goes to zero in the classical limit. 

The inclusion of quantum corrections for only the ideal part of the chemical potential might seem arbitrary, but it 
is thermodynamically consistent, since, Afxpu — >■ as p goes to zero. 

Note that at large Nc such an exact formula is necessary because the simpler formulae, such as the widely used 
[H,!!^ low temperature formula in terms of the Fermi energy e/ = (37r^p/2)^/'^/(2ms) ^ 



(16) 



is not appropriate. If one uses this approximation, the second term of /z, will tend to have an — T^/(jVr.ef) 
contribution which forces fiq away from Aqcd, in contrast to all other effective phase diagrams [35|, |39|, |4J] . The 
problem with this approximation is that it relies on T ^ e/, which, at Nc — c», is not appropriate for any T no 
matter how low. This e/ — >■ scaling is in line with the argument |3l| that at large Nc baryon is a classical object for 
which quantum statistics is inappropriate, and furthermore suggests that quantum effects are irrelevant for collective 
baryon states. Taking into account the fact that heavy atoms tend to make a crystalline structure at low temperature 
(where each atom becomes effectively classical as it is trapped in its location on the crystal), we suggest ( as in 
Hi!, Si, Hi]) an analo gous nuclear ground state in the large Nc limit, as it would provide a physical explanation as to 
why large fermions in the large Nc limit behave like classical objects up to zero temperature. In our world, of course, 
the nuclear liquid is very different from a crystal. 

The resulting quark chemical potential, including all terms of Eq. [T^] and the Nc scaling of a given by Eq. [5] and 
(3 ^ Nc can be seen in the solid lines of Fig. [5l without (panel (a)) and with (panel (b)) the AjipD correction. 
Qualitatively, the diagram looks somewhat different from what we expect the nuclear liquid-gas phase transition to 
look like, but given the roughness of the model presented here this is not so surprising. 

The fact that this diagram looks very similar to the one obtained in [13, is encouraging, through one should not 
forget that the conserved charges examined in [47l l48j are very different from the baryonic charge. The approximately 
linear shape of the diagram should become more familiar with the inclusion of the 7 > curvature term, which would 
also force the chemical potential to smoothly go to a thermodynamically consistent value at T = (currently, Eq. [T^] 
stops being well-defined at a lower temperature since, without a 7-term, P{T > 0) = 0). 

Note also that, generally ps > fns at low temperatures, since the excluded volume plays a much larger effect than 
in our world. Higher order corrections could again fix this, although it is ultimately due to fact that p at the phase 
transition in our world is ^ Aq^^, which, as shown in section UTJ is not natural within large Nc but understandable 
at low Nc- The inclusion oi ppD helps bring the phase diagram quantitatively closer to more realistic calculations 
such as [45|, both at high and low p, but does not change the qualitative structure of the phase transition line. 
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FIG. 5. (color online)The phase diagram in the T — fiq plane for A'^c = 3, 5, 8, 10, 30, 100, without ( panel (a)) and with (panel 
(b)) the AfiFD correction. As Nc increases, the critical points move to the asymptotic limit Nc — oo, represented by the vertical 
dashed line. Both T and fiq are expressed in units of A.qcd- The upper ends of the solid lines correspond to the physical 
critical points. The lower ends of the solid lines correspond to regions where the 7 and higher terms in the VdW expansion 
become dominant (Pressure becomes negative unless 7 is added), and hence is unphysical. 



T 




FIG. 6. (color online) The phase diagram in the T — Hq plane for Nc — 3, 5, 8, 10, 30, assuming P ^ Nc- 

The important conclusion to be drown from Fig. [5] is that, while the Nc 00 curves are qualitatively different 
from realistic nuclear matter and depend on the inclusion of 7 and higher-order terms, the asymptotic Nc — > 00 
limit is independent of these details, and describes the expected shape of the quarkyonic phase transition: The phase 
transition line becomes a vertical line in the fiq — T plane centered around fiq ~ ^qcd- The reason is simply that 
the baryo-chemical potential /zb has a mass component ~ N^, and binding energy components ~ Nj^. When 1^ = 0, 
the latter become negligible, and the chemical potential reaches the limit iiq ~ Aqcd + O (A^^^) at any temperature. 
When v > 1, as we saw in the previous subsection, the phase diagram does not have a physically plausible large Nc 
limit. 

Fig. ini shows what would happen to the fJ- — T phase diagram if the nuclear force term /3 ~ Nc- As can be seen, in 
this case the curvature stays constant at Nc — >■ 00. This is far from surprising, since in this case the binding energy, 
and hence the curvature term on the phase diagram {f{p',T) in Eq ll2p would acquire an iVc-linear component. Fig. 
[6lhowever, also confirms, as the previous sections have shown, that this limit is unphysical at large iVc, since both T 
and n go well above Aqcd soon after Nc — i 



V. DISCUSSION AND CONCLUSIONS 

What are we to make of these findings? If its really true that /3 ~ Nc, as suggested in [i^, we find that the 
Quarkyonic phase and the liquid-gas phase coincide in the large Nc limit. 

The question of whether they also coincide a.t Nc 00 (ie, if we already "discovered the quarkyonic phase" in 
lower-energy nuclear collisions j23l - l30| ). or if the liquid- gas transition is distinct from the Quarkyonic phase in this 
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regime, immediately arises. 

The fact that only one phase transition, with a structure very similar to the one examined here, arises in the 
thermodynamics of charged black holes in AdS spaces [13, 13 ^ thought to be dual to the thermodynamics of strongly 
coupled Gauge theories of large Nc, makes the identification of the quarkyonic transition with the liquid-gas transition 
made in [s^ natural. One should remember, however, that the charges examined in [13,13 (and even in are 
quite distinct from the baryonic charge in QCD. 

Conceptually, the not-so phenomenologically exciting scenario that the liquid-gas transition is all there is is plau- 
sible from universality arguments: The quarkyonic transition arises from a saturation of entropy by inter-baryonic 
interactions [l^, inter-baryonic interactions are in the Van Der Waals universality class, ergo, the quarkyonic phase 
is the liquid gas. 

Furthermore, the question of why does the nuclear liquid look so different from the quarkyonic phase can be sensibly 
answered by looking at the scaling of a: The Van Der Waals terms are universal because they arise when a field giving 
rise to a nucleon-nucleon interaction is integrated out into a nucleon-nucleon interaction term. Since this field, in a 
hot medium, carries entropy, integrating it out under-estimates the entropy of the system. We can now look at Fig. [5] 
to get a sense of how relevant is this undercounting as Nc is varied. If Nc <^ k^d), as in our world, the Pauli exclusion 
principle keeps the nuclear excluded volume at a value significantly larger than Aq^^. In this case, confinement 
suppresses the exchange of colored degrees of freedom between the nuclei, so the entropy carried by whatever forces 
carry nuclear exchange ^ N^. It would also mean that the percolating phase transition studied in [sil. [s^ coincides 
with deconfinement. 

In a Nc 3> k{d) world, however, nuclei touch each other, and colored degrees of freedom can freely percolate 
between them. The entropy carried by these percolating degrees of freedom ~ Nc, and in the large color limit ends up 
overwhelming the total entropy of the system, in much the same way that the electron gas carries most of the entropy 
of a metal (Note that the equilibrium entropy of colored objects ~ Nc even if interaction cross-section between these 
objects is TV^-suppressed. The times cale of equilibration gets longer, but the equilibrium entropy stays the same). In 
this limit, the percolation transition [Sll [52j does not represent deconfinement but the quarkyonic transition, and the 
two, in chemical potential, are separated by A/i ^ NcAqcd- If this scenario is correct, it does indeed imply that 
looking for the Quarkyonic phase in experiments is pointless. 

While the scenario described here is, at the moment, a speculation, this work shows that we do not live in an 
Nc — oo world as far as baryonic matter is concerned, since in our world the number of neighbors dominates over 
the number of valence quarks. The transition from a "low Nc limit" and the "high Nc limit" should happen at 
Nc ^ N]\j ~ 10, and theories with Nc below this scale, such as physical QCD, could look qualitatively different from 
large the Nc calculations. 

The physical manifestations of this transition could conceivably be investigated by varying Nc on the lattice at 
finite chemical potential. Large Nc studies have already been conducted at zero chemical potential [H, [lOl, where 
it was found the thermodynamics is very weakly dependent of Nc already from Nc = 3. Our results suggest that, 
if these studies were extended to finite chemical potential, the results will be very different: At Nc ~ 10 there will 
be a transition from a thermodynamics much like that of our world, with an order-of-magnitude hierarchy between 
nuclear matter and tightly packed baryonic (quarkyonic?) matter, to the "large Nc limit where the two coincide and 
the liquid-gas critical point temperature and chemical potentials J^qcd, and the entropy of the nuclear liquid ~ Nc. 

In the strong coupling expansion (6l| , such a study is already possible with the computational technology available 
today. In the weak coupling lattice limit, the sign problem precludes an exact calculation at finite chemical potential. 
Currently available algorithms i62ri64] become very expensive unless fJ-q/T <^ 1. Together with the added numerical 
cost associated with increasing Nc, this means a true test of A^c-convergence of lattice QCD at finite baryochemical 
potential is still a few years away. 

This difficulty is, however, computational rather than fundamental: Since at large Nc the mixed phase goes from 
very low to very high densities, the lack of convergence with Nc of the critical parameters at finite n should be manifest 
even at relatively small chemical potentials fiq/T <C 1 where the approximation techniques described in [63 - [6^ can 
be used. The computational challenge might be more manageable for 2d lattices, where the "large" number of colors 
is also lower: For the ID QCD ('t'Hooft model [6^]) it will be just Nc = 3 (the fact that any A"c > 2 is in the infinite 
limit in this model might account for its relatively trivial phase diagram [66|), while for 2D QCD it will be Nc ~ 6. 

In conclusion, we have discussed the VdW liquid-gas transition in the context of the large Nc scaling of nuclear 
forces and parameters. We have determined that the VdW model is unsuitable for describing the large Nc nuclear 
liquid-gas transition if the nuclear forces scale as Nc, as generally accepted after [3l|, but will be suitable if these 
forces scale as Af°, as proposed in 44]. In the latter case,intriguingly, the large Nc liquid-gas phase transition coincides 
with the recently proposed Quarkyonic transition [19], although the two are well-distinct at Nc = 3. This hierarchy is 
due to the interplay between two relevant scales, the number of colors Nc and the number of nuclear neighbors Ajv, 
with the second scale larger than the first in the real world but smaller in the large Nc world. A lot of further work, 
both theoretical (lattice simulations, effective theory models, etc) and experimental (low energy scans) is required to 
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interpret these results in a context closer to fundamental QCD. 
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